Quantum cloning of orthogonal qubits 
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An optimal universal cloning transformation is derived that produces M copies of an unknown 
qubit from a pair of orthogonal qubits. For M > 6, the corresponding cloning fidelity is higher than 
that of the optimal copying of a pair of identical qubits. It is shown that this cloning transformation 
can be implemented probabilistically via parametric down-conversion by feeding the signal and idler 
modes of a nonlinear crystal with orthogonally polarized photons. 
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In contrast to classical information, quantum informa- 
tion cannot be copied. This so-called no- cloning theo- 
rem [Q , which is a direct consequence of the linearity of 
quantum theory, states that it is impossible to prepare 
several exact copies (or clones) of an unknown quantum 
state \ip). Although exact cloning is forbidden, one can 
design various quantum cloning machines which produce 
approximate clones. In particular, much attention has 
been devoted to the optimal universal cloning machines 
for qubits, which prepare M identical approximate clones 
out of N replicas of an unknown qubit, and such that the 
fidelity of the clones is state-independent pi. Cloning 
machines for states in a d-dimensional Hilbert space (qu- 
dits) were also investigated p), as well as continuous- 
variable cloning machines for coherent states 

In the limit of an infinite number of clones, the opti- 
mal cloning reduces to the optimal quantum measure- 
ment. In this context, a very interesting observation 
has been made by Gisin and Popescu || who noted that 
the information about a direction in space is better en- 
coded into two orthogonal qubits than in two identi- 
cal ones. If we possess a two-qubit state \ip, ip±) with 
(V'IV'-l) = 5 then we can estimate \ip) with a fidelity 
F x = (1 + l/V3)/2 « 0.789 |, §. This slightly ex- 
ceeds the fidelity of the optimal measurement on a qubit 
pair \ijj,ip), Fu = 3/4. A similar situation occurs for 
continuous quantum variables. Suppose we want to en- 
code a (randomly chosen) position in phase space. A 
possible strategy would be to prepare a pair of coher- 
ent states | a, a), where the real and imaginary parts of 
the complex number a represent the phase-space coor- 
dinates. However, it is actually better to supply a state 
\a,a*} from which a can be inferred via optimal mea- 
surement with a lower error Q. It can also be shown 
that the state \a, a*) gives an advantage when cloning 
coherent states: M identical approximate clones of a co- 
herent state | a) can be prepared with a higher fidelity 
from the state \a, a*) than from |a,cv) 

Motivated by this result, we were led to ask whether 
a similar situation might also occur for qubits. Can M 
clones of qubit be produced from an orthogonal qubit 
pair \ip,ipj_) with a higher fidelity than from an identical 
pair \yj,ip}7 In this Letter, we answer this question by 
an affirmative. We present a universal cloning machine 
acting on an orthogonal qubit pair that approximately 



implements the transformation \ip)\ip±) — > \\j})® M with 
the optimal fidelity. Then, we show that this cloning 
transformation can be implemented probabilistically in 
quantum optics by use of parametric down-conversion. 
Our proposed setup extends the scheme of Simon et al. 
^ by feeding both the signal and idler modes of a non- 
linear crystal with \yj) and |^>±), respectively. 

Let us first provide a simple argument on why we can 
expect the state \ip,il>x) to be better cloned than \ij),vj). 
If we perform an optimal measurement of \tp, ip±), we can 
prepare M identical clones of each with a fidelity F±. 
In contrast, the fidelity of the optimal universal cloning 
machine that prepares M clones from a two-qubit state 
V) is given by F l{ (M) = (3M + 2) /{AM) @. Clearly, 
-Fji(M) < F± for sufficiently large M. Hence, this (non- 
optimal) measurement-based cloning of \vj,vj±) is better 
than the standard cloning of \jp,tp}. 

Let us now seek for a unitary transformation which 
optimally approximates the transformation |V')|V'x) ~~ * 
|V?)® M . Since the set of all states of the form \ij})\ij}±} 
span the whole Hilbert space of two qubits, the most 
general transformation is of the form: 
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where \R) and \Rijk) respectively denote the initial and 
final states of the ancilla, while \M, k) (k = 0, . . . , M) 
denotes a symmetric M -qubit state with k qubits in state 
|0) and M — k qubits in state The arbitrary state of a 
qubit \vj) can be conveniently written as \vb) — d(Q)\0) — 
J2i dio(0,)\i), where the matrix d(£l) is given by 
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with d and <j) denoting the usual polar and azimuthal an- 
gles pointing in direction £1. The linearity of (Q) implies 
that an arbitrary pair of orthogonal qubits transforms 
according to 

W^±) l^outWO) = ^d m (a)d jl {Sl)\M,k)\R iik ). 

ijk 

(3) 

We will measure the quality of the transformation by the 
average single-clone fidelity F±(M). Denoting by Tri )anc 
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the partial trace over the ancilla and all the clones but 
the first one, we get 

F±(M) = /'dO(^|Tr lianc [|'* out (^))(*out(V')l]lV') 

= EE(%''i R ^)4" ( 4 ) 

i' j' k' ijk 

where 

At 3 ( k ' = Y,( n '\Tn[\M,k)(M,k'\]\n) 

(5) 

The coefficients A\^ k k can be considered as matrix el- 
ements of an operator A acting on the space 7i ® K, 
where Ti denotes the Hilbert space of the two input qubits 

and JC denotes the Hilbert space of symmetric states 

•/ -tit 

of M output qubits. Similarly, xljk = (Rijk\Ri> 3 'k') 
define matrix elements of an operator x a l so acting on 
H (8 /C. The formula ([|) for the fidelity thus simpli- 
fies to F±(M) = Tt-h,k[xA]- The operator x uniquely 
represents the completely positive cloning map, which 
transforms operators supported on Ti onto operators sup- 
ported on JC. By definition, the operators A and x are 
Hermitian and positive semidcfinitc, A > and x > 0. 

Of course, the transformation (|l|) should be unitary, 
which reads ^2 k (Ri'j' k I Rijk) =fii>i&j'j- This is equivalent 
to Tr^[x] = lw, where 1^ is the identity operator on 
H. Thus, introducing a set of Lagrange multipliers \\j 
for these unitarity constraints, our problem amounts to 
extremize the quantity W — Ti^ )c[(A — A)%] under the 
constraint x > 0, where A = A (£> TL/c and A is the matrix 
of Lagrange multipliers (Ix; is the identity operator on 
K). Varying W with respect to the eigenstates of the 
operator x-, we get the extremal equation 

[A - A) X = (6) 

for the optimal x- Following |t0]], this equation can be 
further transformed into a form suitable for numerical 



calculation via repeated applications of 

x = A- 1 A x AA-\ \ = {Tv K [AxA]) 1 ' 2 . (7) 

Note that the matrix A > is determined from the uni- 
tarity constraints. 

By numerically solving Eq. (R) for M = 2, ... 15, we 
have been able to conjecture the general analytical form 
of the optimal transformation: 

M 

otj,M\{M - jip±) <g> \(M - j)i/}±, jip), 

3=0 

(8) 

where 

< iy[ 1 , V3(Af-2j) 
J ' [v/2(M+l) v/2M(A/+l)(A/ + 2)J 

(9) 

with | (Af — j)ip±) denoting a totally symmetric state 
of M qubits where j qubits are in state \ip) and M — 
j qubits are in state \ip±). The first M output qubits 
contain the clones of state \ip) while the other M qubits 
contain the clones of ji/'j.) (or anticlones). 

We stress here that the cloning transformation (0) is 
unitary. Since this is by no means obvious from (|8|), let 
us present a proof of this. We can expand any state 
(M — j)ipx) in the basis \M, k) as 

M 

(M - j)j>±) = J2 e i<J - k) +D%(0)\M, k). (10) 

We will not need an explicit expression for the functions 
.Djy (i?) here, but will only use some of their properties. 

Since the functions Djy (i?) are elements of a (real) uni- 
tary matrix, they satisfy the orthogonality relation, 

M 

j2Dllj(mijW=hi. ai) 
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We will also use the following recurrence formula O] , 



(2j-M)Djf j (-d) = {2k - M) costi Djfj(#) + sini? v / (fc + 1)(M - fc) + sind^ k(M - k + 1) D%L ltj (&). (12) 
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For the purposes of the proof it is convenient to apply 
a unitary transformation Uq on the last M qubits at 
the output of the doner and get the state |$ O ut(V0) = 
I® M O Uf M \^ ont (il))). The unitary transformation U Q 
flips the states |0) and |1), |0) -> |1) and |1) -► -|0). 
Thus \(M - j)ip±, jtp) -> (-lp'KM - j)i>*,ji>*i) where 



1^*) = Si d*o I *) - Next we expand |(M - j)i>,jtp±.) an d 
| (M — j)ip* , jipj_) in the basis |M, fc) according to Eq. 
(fiOj) , and then utilize the recurrence formula (|12|). Fi- 
nally, we can carry out the sum over j with the help of 
Eq. (pi]), resulting in 
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where the coefficients au and &m read 
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The four terms on the right-hand side of Eq. (|13|) are 
proportional to the output states for the four input ba- 
sis states 1 1} , 1 10} , 1 00) , and |11), respectively. Con- 
sequently, it is easy to prove that the transformation 
\ip,ip±) — ► |3>out(V0) preserves scalar products, hence is 
unitary. 

Let us now calculate the fidelity of the clones. We 
can see from Eq. (||) that the cloning machine preserves 
the symmetry of the input state \i/),ipx), so the clones 
of both states \tp) and \ip±) have the same fidelity. This 
state-independent single-qubit fidelity can be obtained 
by summing a series, 
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After some algebra, we arrive at the expression 
1 
2 
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(14) 



(15) 



We are now able to compare this fidelity to that of the 
optimal doner for a pair of identical qubits Fj | (M) : for 
M < 6, we have F\\{M) > F±(M), while Fj|(6) = F±(6) 
and F±(M) > F\\(M) for M > 6. Thus, the doner (§) 
outperforms the standard doner for M > 6. We note 
also that for M — > oo, the fidelity F±(M) tends to the 
optimal measurement fidelity F±, as expected. 

The optimality of the doner can be proved with the 
help of techniques adapted from the theory of semidcf- 
inite programming |L2|. We observe that the trace of 
Lagrange multiplier A provides an upper bound on the 
achievable fidelity. If A Cg> Ik — A > then it holds 
for any X that Tr n ,ic[x* ® Ik] > Trn,lc[x A ]- lt fol " 
lows from the unitarity constraint Tr^[x] = 1-h that 
T^H,lc[x^ ® Ijc] = Tr«[A] does not depend on \- Thus 
it holds that Tr^fA] > Tr^x; • From the numerical 
solution of Eqs. (@) we have 'in basis |00), |11), |01), 1 10) , 



A = 



F ± (M) 
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(16) 



The block-diagonal matrix A® t/c — A is positive semidef- 
inite and has three different eigenvalues which read fix = 



fc=0 



(13) 



M+2 
3M ' 



and fis = 0. Since the up- 



J_ /M+2 _ 1 

12 V 3M ' ^ 2 — 3 

per bound Tr^[A] = is saturated by our cloning 

machine, we conclude that our doner is optimal. 

In the rest of this paper, we will show that the cloning 
transformation (^J) can be implemented probabilistically 
via stimulated parametric down-conversion. The exper- 
imental setup under consideration is shown in Fig. 1. 
This scheme is a straightforward extension of the setup 
suggested by Simon et al. Q where the qubits are rep- 
resented by the polarization state of photons. We can 
identify |0) with vertical polarization and |1) with hor- 
izontal polarization states. In optical parametric down- 
conversion, a 'blue' photon can split into a pair of 'red' 
photons. Traditionally, these daughter photons are re- 
ferred to as signal and idler, respectively. In our setup, 
three nonlinear crystals Ci, C2, C3 are pumped by a 
strong laser beam. In crystals C\ and C2, pairs of pho- 
tons can be produced, so we can verify the presence of 
signal photons by detecting the idler photons emerging 
from C\ and C2. If a single idler photon is detected on 
each side, then we have one signal photon in each beam. 
The states of these two photons can be manipulated with 
the help of phase shifters and polarization rotators in or- 
der to prepare the desired input state \ip,ip±). The two 
photons then feed the signal and idler modes of a third 
nonlinear crystal C3, where M clones are generated due 
to the stimulated parametric down-conversion. 
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FIG. 1: Setup for the cloning of orthogonal qubits via stimu- 
lated parametric down-conversion. For a detailed description, 
see text. 

In the limit of strong coherent pumping, the effective 
Hamiltonian describing the interaction in C3 can be writ- 
ten as follows P|, 



H = iTig(a vl a\ 



H2 



l Hl u 'V2! 



h.C. 



(17) 
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where a* vi and a' m denote bosonic creation operators for 
photons in the first mode with vertical (V) or horizon- 
tal (H) polarization, and similarly a|, 2 and a} H2 are cre- 
ation operators for photons in the second spatial mode. 
The constant g denotes the parametric gain. The time 
evolution is thus governed by the unitary transformation 
U = ex-p(—iHt/K). With the help of the disentangling 
theorem, we can write the operator U in a factorized form 

{J = e ra vi a L2( cos l 1 7) _ ( a vi av ' 1+a H2 a «2 + l) e -ravia H 2 

xe~ ra Hi a ^2(cosh7)~( a Hi aH1+a ^2 a ^2+i) e raHiai^ 



where 7 = gt and T = tanh.7. The Hamiltonian (|1^ 
has the important property of being invariant under gen- 
eral simultaneous SU(2) transformations on the polar- 
ization vectors (ay, an) for modes 1 and 2 |J. It is 
thus sufficient to consider the evolution of a basis state 
|l)vi|0)iji|0)y2|l)ff2 ( a single vertically polarized pho- 
ton in mode 1 and a single horizontally polarized photon 
in mode 2) which represents the input state ifix) = 
1 01) . Making use of the factorized form of U , we obtain 
the state at the output of the crystal C3 in the form 

00 M 

£ r^(i - r 2 ) [(m - - r 2 ) - r 2 ] 

M=0 j=0 

x \M - j)vi \j)m \j)v2 \M - j) H2 , (18) 

where \k)i with I — VI, HI, V2, H2 denote the usual 
Fock states. For a fixed number M of photons in each 
mode 1 and 2, the output state (|l^) closely resembles the 
output state of the universal cloning machine m) with the 
coefficients a jjM (T) w [{M - - T 2 ) - T 2 ] (-l)J. If 
we measure the number of photons in mode 2 and detect 
M photons, then we know that M photons representing 
M approximate clones of the input qubit \ip) are present 
in mode 1. In order to calculate the fidelity of these 
clones, we insert the properly normalized o^m^) into 
formula (14). After some algebra, we obtain 

x 3y 2 - 2y(2M + 1) + §M(M + 1) , s 
F (M, y) = , ^ (19) 



where we have introduced y — T 2 /(l — T 2 ) = sinh 2 7 
for notational convenience. The cloning fidelity thus de- 
pends on the parametric gain 7, so we must optimize 
this gain in order to achieve the highest possible fidelity. 
Upon solving dF ^ I ' y ^ ! — for y, we find that 



2/opt 



M 
T 



M(M + 2) 



(20) 



By inserting y opt into Eq. (pL 9|) , we recover the optimal 
fidelity (|l5|). Furthermore, it can be verified by direct 
calculation that with the optimal gain, the postselected 
M -photon state at the output of the crystal C3 coincides 
with the output of the cloning machine (||) . 

This approach of cloning based on down-conversion can 
be further extended to the approximate realization of the 
general cloning transformation \^)® N \iP±)® N ' -> \if;)® M . 
For N' = 1, we have been able to derive the optimal 
fidelity for any N and M > N by a similar calculation, 



s N+l 3(N - 1) + JP/(N + 2) , x 
F±(N,M) = —^— + ^ n '^ T ■ ;\ '- (21) 



N + 3 



2M(N + 3) 



6y 2 - 6My + M(2M + 1) 



withP= (N-l){N 2 -15N-18)+8M(N+l){M+3-N). 
It can be checked that there is again a value of M above 
which this doner outperforms the standard (iV+1) — * M 
doner. For large N, however, the advantage becomes 
marginal. 

In summary, we have designed a universal cloning ma- 
chine for orthogonal qubit pairs, and have shown that it 
achieves a higher fidelity for M > 6 than the standard 
doner for identical qubits. We conclude that the advan- 
tage of orthogonal qubits over identical qubits that was 
discovered in the context of measurement also extends to 
cloning. 
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